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Abstract. - We investigate the momentum distribution function near the Mott-Hubbard 
transition in the one-dimensional ti — £2 Hubbard model (the zig-zag Hubbard chain), with 
the density-matrix renormalization-group technique. We show that for strong interactions 
the Mott-Hubbard transition occurs between the metallic-phase and an insulating dimerized 
phase with incommensurate spin excitations, suggesting a decoupling of magnetic and charge 
excitations not present in weak coupling. We illustrate the signatures for the Mott-Hubbard 
transition and the commensurate-incommensurate transition in the insulating spin-gaped state 
in their respective ground-state momentum distribution functions. 



Introduction - Many aspects of the low-energy physics of an electronic system are influ- 
enced by the shape of its Fermi surface and the occupation of nearby states. Finite tempera- 
ture, disorder or strong correlations may change the Fermi-surface geometry and or topology 
and induce magnetic and other instabilities. 

The investigation of the interplay of magnetic interactions, in particular the existence 
of incommensurate phases, with charge excitations emerges as one of the central issues in 
the physics of low-dimensional electronic systems [1-7]. The nature and mechanism of the 
Mott-Hubbard transition, as one of its most striking manifestations, has been subject to 
intense investigation for many years. It is therefore important to study the renormalization 
of individual Fermi-surface sections under the influence of electron-electron correlations in 
competition with frustrating interactions. 

Recent investigations suggested a spontaneous, interaction induced deformation of the 
Fermi surface of the 2D (extended) Hubbard model close to half filling, indicating in part a 
violation of the Luttinger sum-rule [8-10]. The possibility of a relevant Fermi-surface renor- 
malization near the Mott-Hubbard (MH) transition in one-dimensional generalized Hubbard 
models has been raised on the basis of a DMRG-study [11], drawing on arguments from results 
obtained for the Fermi-surface flow by RG [12,13]. 

In this letter we examine this Mott-Hubbard transition in one of the prototypical frustrated 
one-dimensional models: the half-filled zig-zag Hubbard ladder. For appropriate choices of the 
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parameters this model describes either the low-energy properties of Hubbard-ladders [12] or 
half- filled edge-sharing double-chain materials like SrCuC>2 [14] or L1V2O5 [15] for which the 
next-nearest neighbor hopping is expected to be substantial. The study of this system per- 
mits both the consideration of incommensurate phases and strong interactions. The former 
are difficult to study in the limit of infinite dimensions, where part of our present under- 
standing of the Mott-Hubbard transition orginates. The latter can now be treated adequately 
for one-dimensional systems using the density-matrix renormalization group (DMRG), which 
emerged in the last decade as a reliable tool to investigate the electronic structure of quasi 
one-dimensional systems [16,17]. 

We find no renormalization to perfect-nesting of the Fermi-surface in the metallic state 
as the Mott-Hubbard transition is approached. We explicitly demonstrate that the dimer- 
ized state on the insulating side of the Mott-Hubbard transition has gaped incommensu- 
rate spin-excitations, which only later give way to an insulating phase with commensurate 
spin-excitations. Our results indicate that perfect-nesting is not a prerequisite for the Mott- 
Hubbard transition at finite values of the interaction. The opening of the charge gap decouples 
from changes in the nature of magnetic excitations, in stark contrast to the weak-coupling 
scenario, where Umklapp-scattering becomes relevant at perfect nesting and leads to the open- 
ing of a charge-gap. We report the characteristic signatures of the incommensurate insulating 
state in momentum-distribution function for future experimental analysis. 

Model - The Hamiltonian of the zig-zag Hubbard ladder is given as: 

H = Yl tAn (4+An,<T C n, CT + hx -) + U ^2 n t n h (!) 

A^l,2 

where the a n IT (c na ) are Fermion-creation (destruction) operators on site n and spin a =f, J, 
and n a = 4,<r c n )0 ~ 

Let us first discuss a few known properties of the phase-diagram at half filling, compare 
Fig. ^ We follow Balents and Fisher [18] and denote with 

CnSm 

a phase with n/m gapless charge/spin modes. A Mott-Hubbard transition of type C0S1- 
C2S2 is predicted [12] by RG at half-filling for t 2 = ti/2. This prediction should be valid at 
infinitesimal U jt\ . 

In the limit of large U the model (JIJ transforms to the J\ — J2 chain 

Hj = Jati S„ • S n+ An j (2) 

Ati = 1,2 

with Ja-,i = it\ n /U. The J\ — J2 model spontaneously dimerizes [19] for big enough J2. The 
critical value for a = J2/ Ji can be determined by examination of the nature of the lowest-lying 
excitation for finite clusters [20]. In the spin-fluid state C0S1 it is a triplet, in the dimerized 
state C0S0 it is a singlet (the dimerized state is doubly degenerate in the thermodynamic 
limit). One finds [20]: a c = 0.2411 = (0.491) 2 . 

Torio et al. [21] have extended the level-crossing study used to determine a c for the J\ — J2 
model to finite U jt\ see the dashed line in Fig. ^ They propose that the dimerization line 
C0S1-C0S0 extends until U/t\ — > and connects to t 2 = h/2. This would consequently 
invalidate the weak-coupling RG-prediction of a Mott-Hubbard transition of type C0S1-C2S2. 
It seems, presently, more likely that the transition is of COSO-CnSm type, with n > 0. 
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Fig. 1 - Phase diagram for the half-filled t\ — ti Hubbard model as explained in the text. Included 
are the DMRG-result for the Mott-Hubbard transition at U = 4£i (filled diamonds) and for the 
commensurate-incommensurate transition in the insulating state (open diamonds). The dashed- 
dotted lines are guides to the eye. The states CnSm predicted by weak-coupling RG are denoted at 
the bottom of the phase diagram. 



The Majumbdar-Gosh point J2 = Ji/2, also called disorder-point, has an exact valence- 
bond dimer ground-state. For a > 0.5 = (0.707) 2 the short range spin-spin correlations 
become incommensurable. Due to the absence of long-range order [22], the peak in the 
static structure factor S(q) moves away from q = 7r only at the Lifschitz-point, which can be 
determined via DMRG [23] as a L = 0.52063 = (0.7215) 2 . 

The classical spin- wave solution to (J2J yields a long-ranged spiral ground-state for a > 0.5 
[24,25]. A long-ranged ordered states is unstable towards quantum-fluctuation in ID, they 
become short-ranged. For the J\ — J2 model the quantum fluctuations lead to a dimerized 
state with finite dimerization and short ranged spiral (incommensurate) spin-spin correlations. 
They show up as a peak in S(q) for a incommensurate wavevector q and width l/£, where £ 
is the correlations length. There has been no study so far of the extension of the dimerized 
phase with incommensurate excitations, realized for a > 0.5 in the limit of large U jt\ to finite 
values oiU/t\. Here we propose that it connects to the metallic state. 

Also included in Fig. are the predictions for the critical U c for the Mott-Hubbard tran- 
sition by Torio et al. [21] and by Aebischer et al. [26]. It has turned out to be very difficult to 
determine numerically the location of this transition from estimates of the charge gap, extrap- 
olated to the thermodynamic limit, due to the fact that the charge gap is exponentially small 
near the Mott-Hubbard transition. Here we find that a quite accurate lower bound for the 
the Mott-Hubbard transitions can be obtained form the study of the momentum distribution 
function. 

Method - In the last decade the density matrix renormalization group [16] emerged as a 
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Fig. 2 - Illustration of the momentum-distribution function n(k) for L = 48 sites, U = Ati and 
t 2 /ti = 0.6, 0.8. DMRG results for iV T = iVj. = L/2 (filled circles) and for iV T = L/2 + 2, ATj. = L/2- 2 
(triangles-up: majority-spin, triangles-down: minority-spin). The lines are fits by Eq. @J. 



reliable tool for the study of electronically one-dimensional models [17]. For strongly inter- 
acting fermionic systems, the evaluation of n(k) as the Fourier transform of the correlation 
function: 

2 L 

™<r( fc ) = ^ 2l cos(fc(n-n')) (4 iCT c„ /)0 .) (3) 

n,n' — 1 

is numerically difficult and costly in the framework of the DMRG [16,17], in particular for 
periodic boundary conditions [27]. Here we report results from DMRG calculations on half- 
filled chains of length L=48, 80 using up-to 300 DMRG states, which we found sufficiently 
accurate to determine n(k). 

To obtain estimates for the Fermi wavevectors we have analyzed the momentum distribu- 
tion function, see Eq. obtained by DMRG, via two smoothed step functions at &fi and 
fcf?2 with respective widths pi and pi [11]: 

k — kpi k — kp2 

n(k) = ciq + a\ atan h a>2 atan . (4) 

Pi " P2 

The quality of fits by J3J to n{k) is illustrated in Fig. We will use the such obtained 
estimates for the Fermi- wavevectors kpi and kp2 in the further analysis. 

Results - In order to elucidate the properties of the Mott-Hubbard transition [26] and 
of the incommensurate-commensurate transition in one dimensional frustrated systems, we 
studied the momentum-distribution function for a range of parameters that crosses all relevant 
phase transition-lines of the phase diagram of the ID ti — t 2 Hubbard model [28]. In Fig. El 
we present our estimates for kpi as a function of ^Ai for U = 4ti, in comparison with the 
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Fig. 3 - For U — 4ti and L = 48,80 sites the DMRG results for JtFi and the difference k F 2 — k F i 
(both in units of n) as obtained by fitting the DMRG-results for n(fc) by Eq. (|IJ, compare Fig. [5] The 
lines are guides to the eye. Below t^/tx ~ 0.74 two things happen: k F i renormalizes substantially 
with respect to its U = 0- value (dashed line) and k F 2 — k F \ deviated from its Luttinger-sum-rule- value 
of 0.5. 

result for {7 = 0. We note a substantial renormalization of k F \ towards larger values below 
t 2 /h « 0.74. 

In the metallic state the Luttinger sum-rule states that the total volume of both Fermi- 
seas together equals the number of electrons. At half-filling this statement is equivalent with 
fci?2 — kpi = 7t/2. Our DMRG-data for (k F2 — k F i) /ir presented in Fig.[3]mdicates a violation 
of the Luttinger sum-rule below *2 / ^i ~ 0.74. We therefore conclude that the Mott-Hubbard 
transition occurs in the vicinity of this point (as indicated in Fig. ^ by the filled diamond in 
Fig. ^) and that the system is insulating for t%lt\ < 0.74. Due to the exponentially small 
gap in this region our result for the position of the Mott-Hubbard transition is in fact a lower 
bound (in terms of t 2 /t\) for the exact transition point, as we would not be able to resolve 
numerically a possible exponentionally small departure from the Luttinger-liquid sum-rule for 
ta/ti > 0.74. 

Next we turn to an analysis of the nature of the spin excitation in the vicinity of the Mott- 
Hubbard transition. The spin-excitation in the insulating state correspond to renormalized 
particle-hole excitations with a spin-flip. They correspond therefore to states with iVf = 
L/2 + AN/2 and N± = L/2 — AN/2. Using for numerical convenience A7V = 2 we have 
calculated the respective Fermi-wavevectors for the majority /cfij an d minority kpi.\, spins. 
We define the normalized shift for the Fermi-wavevector as 

Ak F1 = {k F1A -k Fld )L/(TrAN) . (5) 

The normalization in Eq. (J5J is chosen such that Ak F i = 1, independent of system size L, 
for case of a single Fermi-sea. In the limit ^Ai — > oo the Hubbard zig-zag chain has two 
equally large Fermi-seas, with equal Fermi-velocities, and Akpi —> 0.5 in this limit, as we 
have verified numerically for large t 2 /ti. 

We have obtained Ak F i by fitting the respective momentum distribution functions for 
the majority and minority spins, see Fig.|2 by Eq. JIJ. The results are presented in Fig. ^ 
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Fig. 4 - The normalized shift A/cfi, see Eq. for L = 48,80 sites and U = 4ti. The lines are 
guides to the eye. For comparison, we have reproduced the results for &f2 — kpi from Fig. 

We notice that Akp\ w 1 for small t^/ti, implying commensurated spin excitations. A well 
defined kink in Akpi at tz/tx w 0.64 — 0.66 indicates a second-order transition to a state with 
incommensurate spin-excitations (denoted by the open diamond in Fig.^). 

These results imply that the Mott-Hubbard transition takes place for Fcrmi-wavevectors 
at arbitrary, incommensurate values, and the Fermi- wavevectors do not renormalize towards 
perfect nesting fcj? i — ■> 7r/2 and &F2 — * 7r. In weak-coupling, Umklapp-scattering becomes 
relevant at perfect nesting and leads to the opening of a charge-gap. Our results indicate, 
that effective perfect-nesting is not necessary for the Mott-Hubbard transition for finite values 
of the interaction-strength. 

Discussion - The interplay between magnetic order and carrier mobility in strongly inter- 
acting systems, with the Mott-Hubbard transition as one important manifestation, has long 
been subject to intense investigation. For the two-dimensional (2D) case, important for exper- 
imental realizations, exact analytical or numerical techniques are still lacking. It is known that 
in 2D the commensurate ordered phase for the nearest-neighbor Hubbard model is unstable 
with respect to helical fluctuations [1,2] either upon doping or in the presence of frustrat- 
ing next-nearest-neighbor interactions. Other, even more exotic phases [3,4] have also been 
proposed, in particular in the context of studies into the mechanism of high-temperature su- 
perconductivity. Our investigation of the one-dimensional case offers new insight into driving 
forces behind the MH transitions: nesting, i.e. the divergence of the magnetic susceptibility 
due to special features of the Fermi surface geometry emerges as sufficient, but apparently not 
a necessary condition of the MH transition [29]. This result is inaccessible by weak coupling 
theory and indicates the existence of a generic mechanism for charge carrier freezing that is 
independent of dominant Umklapp scattering processes. Spin-charge separation in Luttinger 
liquids may enhance this effect in one-dimension, but similar strong-coupling scenarios must 
be present also in higher dimensions to explain the MH transition between incommensurate 
insulators and metals in doped or frustrated 2D systems. 
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